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ABSTRACT 

Using Suzaku and the Rossi X-ray Timing Explorer ( RXTE), we have conducted a series of four simultane- 
ous observations of the galactic black hole candidate Cyg X-l in what were historically faint and spectrally hard 
“low states”. Additionally, all of these observations occurred near superior conjunction with our line of sight 
to the X-ray source passing through the dense phases of the “focused wind” from the mass donating secondary. 
One of our observations was also simultaneous with observations by the Chandra-High Energy Transmission 
Grating (HETG). These latter spectra are crucial for revealing the ionized absorption due to the secondary’s 
focused wind. Such absorption is present and must be accounted for in all four spectra. These simultaneous 
data give an unprecedented view of the 0.8-300 keV spectrum of Cyg X-l, and hence bear upon both corona 
and X-ray emitting jet models of black hole hard states. Three models fit the spectra well: coronae with thermal 
or mixed thermal/non-thermal electron populations, and jets. All three models require a soft component that 
we fit with a low temperature disk spectrum with an inner radius of only a few tens of GM /c 2 . All three models 
also agree that the known spectral break at 10 keV is not solely due to the presence of reflection, but each gives 
a different underlying explanation for the augmentation of this break. Thus whereas all three models require 
that there is a relativistically broadened Fe line, the strength and inner radius of such a line is dependent upon 
the specific model, thus making premature line-based estimates of the black hole spin in the Cyg X-l system. 
We look at the relativistic line in detail, accounting for the narrow Fe emission and ionized absorption detected 
by HETG. Although the specific relativistic parameters of the line are continuum-dependent, none of the broad 
line fits allow for an inner disk radius that is > 40 GM/c 2 . 

Subject headings : accretion, accretion disks - black hole physics - radiation mechanisms :non-thermal - X- 
rays:binaries 


1. INTRODUCTION 

There is currently significant debate as to the physical 
mechanisms responsible for the continuum of X-ray spec- 
trally “hard states” of black hole candidates (BHC) accret- 
ing in binaries. This debate ranges from the broader issue of 
whether or not there is a significant contribution to the X-ray 
band from an outflow or jet, to more narrowly focused issues 
within given classes of models. For instance, the hard X-ray 
emission has traditionally been attributed to a Comptonizing 
thermal corona (Eardley, Lightman & Shaipro 1975; Shapiro, 
Lightman & Eardley 1976; Sunyaev & Trumper 1979; Dove 
et al. 1998, and references therein). These earlier works gen- 
erally favor a scenario where the corona lies central to a trun- 
cated outer thin disk. However, if the corona is driven out- 
wards by radiative pressure (e.g., Beloborodov 1999) could it 
instead overlay the inner disk? Can the optically thick, ge- 
ometrically thin disk extend inward nearly to the innermost 
stable circular orbit (Miller et al. 2006)? Is this disk cold 
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(peak temperatures of a few hundred eV), or can it instead 
be hot (near a keV, i.e., Wilms et al. 2006)? Is the hard 
state corona comprised primarily of electrons with a thermal 
population (Poutanen & Vurm 2009), or can it have a sub- 
stantial contribution from a non-thermal electron population 
(Ibragimov et al. 2005)? Does the bulk motion of the flow 
play a role in Comptonizing the spectrum (Shaposhnikov & 
Titarchuk 2006; Laurent & Titarchuk 2007)? Alternatively 
could the X-rays be comprised of a combination of direct syn- 
chrotron and synchrotron self-Compton (SSC) emission from 
a jet (Markoff, Nowak & Wilms 2005; Maitra et al. 2009)? 
Contributing to the debate, however, is the fact that many of 
the above cited models, especially when considering solely 
Rossi X-ray Timing Explorer (RXTE) data in the 3-200 keV 
range (or even narrower energy ranges), describe the hard 
state spectra nearly equally well. 

To study these issues, over the past decade we have been 
using a series of pointed, approximately bi-weekly, RXTE 
observations of the BHC Cyg X-l performed simultane- 
ously with 15 GHz radio observations by the Ryle telescope. 
Cyg X-l holds much promise for exploring the current range 
of questions listed above owing to its persistently bright X- 
ray flux (including both the hard and soft states, it varies be- 
tween 200-600 mCrab in the 1 .2-1 2 ke V band covered by the 
RXTE -All Sky Monitor) and its correlated radio/X-ray spec- 
tra (see Wilms et al. 2006, and references therein). Extended 
radio emission even has been imaged in Cyg X-l (Stirling 
et al. 2001). This campaign has already provided the spectra 
for some of the Comptonization models (Wilms et al. 2006), 
and jet models (Markoff, Nowak & Wilms 2005) discussed 
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above. Additionally, these data have been used to to study 
the correlation of the X-ray and radio spectral properties on 
both long time scales (Pottschmidt et al. 2000; Pottschmidt 
et al. 2003; Wilms et al. 2006), and short time scales (Gleiss- 
ner et al. 2004, Wilms et al. 2007, Bock et al. 2010, Bock et 
al., in prep.). Thus they comprise a strong data set for ad- 
dressing not only the details of the X-ray spectrum, but also 
the connection to the hard state jets which are known to dom- 
inate the radio through near-infrared emission. Simultaneous 
radio/X-ray flaring has also been detected in Cyg X-l dur- 
ing this extended campaign (Fender et al. 2006, Wilms et al. 
2007), lending support to the hypothesis of X-ray emission by 
the jet. 

RXTE spectral data alone, however, do not allow us to 
break the current existing “theoretical degeneracy” in the ori- 
gin of the X-ray spectrum. The statistically best fits are in 
fact obtained with purely empirical simple broken powerlaws 
with a break occurring between 9-12keV, and an exponen- 
tial cutoff occurring at >20keV, to which a broad, « 6.4keV 
gaussian line is added (Wilms et al. 2006; Nowak et al. 2005). 
The latter component is likely attributable to a relativistically 
broadened Fe Ka line (Reynolds & Nowak 2003, and refer- 
ences therein); however, its parameters are dependent upon 
the assumed continuum model (Wilms et al. 2006). More 
physically motivated Comptonization and outflow-dominated 
models (Markoff, Nowak & Wilms 2005; Wilms et al. 2006) 
can describe the same spectra almost as well. However they 
must introduce additional, albeit plausible, physical compo- 
nents (e.g., relativistic smearing) to recover the simple spec- 
tra of the broken powerlaw description. Thus, there is some 
amount of ambiguity when correlating detailed spectral fea- 
tures vs. continuum properties, e.g., disk reflection vs. coro- 
nal compactness/hardness, as the detailed features systemati- 
cally depend upon the underlying broad-band continuum. 

Despite these problems in finding a truly unique spectral 
model, when considering multiple observations taken over a 
wide range of luminosities and spectral hardnesses, spectral 
correlations arise that are robust and persistent across a vari- 
ety of these theoretical characterizations (Wilms et al. 2006). 
Using the broken powerlaw models as a simple description 
of the X-ray spectra, we have found that when the 2-10keV 
photon index T i < 2.2 there is a positive correlation between 
X-ray and radio flux, whereas for Ti > 2.2 there is a negative 
correlation between the X-ray and radio flux. We therefore 
use the value of T i ps 2.2 as the canonical division between 
the spectrally “hard state” and the spectrally “soft state” (see 
Remillard & McClintock 2006). Additionally, as the spectra 
become harder, exponential cutoffs tend to become less sig- 
nificant (see also Motta, Belloni & Homan 2009). 

Our previous RXTE spectral studies of Cyg X-l have been 
limited in two respects: the low spectral resolution of RXTE 
{E/XE 6 at 6keV), and the inability to measure spectra at 
^ 3 keV. In this work, we turn to a set of four Suzaku ob- 
servations that we performed simultaneously with our RXTE 
campaign to enhance our Cyg X-l studies in several crucial 
ways. First, Suzaku has large effective area at soft X-ray en- 
ergies. In this work, we consider spectra down to 0.8 keV, 
which allows for the possibility of measuring the “seed pho- 
ton spectrum” in Comptonization models (§5.1), or judging 
the relative contribution of synchrotron versus the disk radia- 
tion in jet models (§5.2). Suzaku also has excellent resolution 
in the Fe Ka line region (E / AE « 50), which allows us to 
separate narrow from relativistically broadened line features 
(§4.3, 5.3). Third, Suzaku measures the Cyg X-l hard X- 


TABLE l 

Log of Cyg X-l Observations 


Date 

(yyyy-mm-dd) 

Spacecraft/ObsID 

Instrument 

Exposure 

(ksec) 

2006-10-30 

Suzaku/401059010 

XIS 0-3 

35. 0 a 

32.1 



HXD-PIN 

27.7 

24.9 



HXD-GSO 

27.7 

25.8 


RXTE/80110-01-13 

PCA 

8.0 b 

8.0 



HEXTE-A 

3.1 

3.1 



HEXTE-B 

2.5 

2.5 

2007-04-30 

Suzaku/40207201 0 

XIS 0,1,3 

34. 0 a 

22.4 



HXD-PIN 

40.2 

27.6 



HXD-GSO 

40.2 

27.1 


RXTE/92090-01-16 

PCA 

14.6 b 

1.5 



HEXTE-A 

4.6 c,d 

0.3 



HEXTE-B 

3.0 

0.4 

2007-05-17 

Suzaku/402072020 

XIS 0,1,3 

22. 3 a 

12.3 



HXD-PIN 

32.6 

10.9 



HXD-GSO 

32.6 

10.9 


RXTE/ 92090-01-17 

PCA 

2.0 b 

0.0 



HEXTE-A 

1.9°’ d 

0.0 



HEXTE-B 

0.9 

0.0 

2008-04-19 

Suzaku/40306501 0 

XIS 1,3 

16. 9 a 

7.2 



XIS0 

34. 0 e 

0.0 



HXD-PIN 

29.0 

10.1 



HXD-GSO 

29.0 

10.1 


RXTE/93 1 20-01 -01 

PCA 

21. 5 b 

7.2 



HEXTE-A 

17.9° 

7.1 



HEXTE-B 

10.8 

4.6 


Chandra/%525 

HETG 

29.4 f 

11. 1 8 


Note. — Exposure times are after initial good time filtering (left), and 
after color/intensity time filtering (right). 
a Summed exposure times for all listed XIS detectors. 
b Summed exposure intervals, not weighted by the fraction of operating Pro- 
portional Counter Units (PCU). 
c HEXTE-A cluster in fixed position (no rocking). 
d Evidence for HEXTE-A cluster exposure time anomaly. 
e XJS 0 run in continuous readout mode. 
f Summed exposure times for all 1 st order spectra. 

8 Summed exposure times for 1 st order HEG spectra. 

ray spectrum up to « 300 keV (§4.1), providing further con- 
straints on Comptonization and jet models. 

The outline of this paper is as follows. In §2, we describe 
our full set of observations and our data reduction procedures. 
Due to maintenance and upgrade of the Ryle radio telescope 
during the construction of the Arcminute Microkelvin Imager 
(AMI), no simultaneous radio measurements are available for 
these Suzaku/RXTE observations. Instead, we use previous 
observations to estimate the radio fluxes (§2.4). For one of our 
observations simultaneous Chandra-HETG data are available 
(§2.3). These Chandra data become crucial in all of our anal- 
yses as they help elucidate the spectral variability associated 
with the observed lightcurve behavior, as discussed in §3. We 
present simple phenomenological descriptions of the spectra 
in §4, including a description of the composite line profile 
(§4.3). Comptonization and jet models are presented in §5, 
along with further discussions of the implied relativistic lines 
(§5.3). We summarize our findings in §6. 

2. OBSERVATIONS AND DATA ANALYSIS 

Prior to 2009 April, there have been five Suzaku observa- 
tions of Cyg X-l. The four most recent of these observations 
are discussed in this work. The first Suzaku Cyg X-l ob- 
servation occurred in 2005 October (Makishima et al. 2008). 
During this observation Cyg X-l was in a relatively bright 
hard state (the RXTE -ASM count rate was « 30cps), the 
spacecraft aimpoint was placed on the Suzaku -X- ray Imag- 
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ing Spectrometer (XIS) detectors, and Suzaku was run in a 
data-mode with 1 s integrations per CCD exposure frame. As 
discussed by Makishima et al. (2008, see also the Appendix), 
this long exposure led to both telemetry dropouts and severe 
photon pileup on the detectors. As these issues require a more 
complex analysis (Makishima et al. 2008), we do not consider 
these data further. 

The four Suzaku observations discussed here occurred dur- 
ing times when the RXTE -ASM flux ranged from 12-23 cps, 
the spacecraft pointing was set to the Hard X-ray Detector 
(HXD; this slightly reduces the flux on the XIS detectors), 
and the CCD exposure frame integration times were set to 
« 0.5 s. Thus, our observations do not suffer from telemetry 
dropouts and they are less severely affected by photon pileup. 
All four of these observations occurred simultaneously with 
RXTE observations. The last of these observations also oc- 
curred simultaneously with observations by every other X- 
ray/soft gamma-ray instrument flying at that time (Hanke et 
al., in prep.). Here we only consider the simultaneous Chan- 
dra-High Energy Transmission Grating. An observing log is 
presented in Table 1 . 

2.1. Suzaku Analysis 

The Suzaku data were reduced with tools from the HEA- 
SOFT v6.8 package and calibration files dated 2009 Septem- 
ber 25. The instruments on Suzaku (Mitsuda et al. 2007) are 
the X-ray Imaging Spectrometer (XIS; Koyama et al. 2007) 
CCD detector covering the « 0.3-1 OkeV band, and the Hard 
X-ray Detector (HXD; Takahashi et al. 2007) comprised of 
the PIN diode detector (PIN) covering the « 10-70keV band 
and the gadolinium silicate crystal detector (GSO) covering 
the « 60-600 keV band. The XIS has four separate detec- 
tors, XIS 0-3, with XIS 1 being a backside illuminated CCD. 
XIS 2 was lost due to a micrometeor hit in late 2006, and thus 
was available only for the first observation (see Table 1). For 
our fourth observation, XIS 0 was run in continuous readout 
mode, which is not yet fully calibrated, so we do not include 
these data in this work. 

In preparing the XIS spectra, we first corrected each detec- 
tor for Charge Transfer Inefficiency (CTI) using the xispi 
tool, and then reprocessed the data with xselect using the 
standard xisrepro selection criteria. Due to thermal flex- 
ing of the spacecraft, the attitude of the Suzaku spacecraft 
exhibits variability over the course of the observations and 
therefore the image of the source is not at a fixed position 
on the CCD. Standard processing reduces this variability and 
improves the PSF image (Uchiyama et al. 2008); however, 
with the standard tools it is not possible yet to correct fully 
the blurring caused by the varying attitude. As described in 
the Appendix, for bright sources a better reconstruction of the 
attitude solution, and thus narrower PSF images, can be ob- 
tained using the aeatt cor . si software. We produced such 
an improved image and then used the pile_estimate . si 
S-Lang script described in the Appendix to estimate the de- 
gree of pileup in the observations. For the spectra described 
in this work, using the time filter criteria described in §3, the 
center of the PSF images could be affected by as much as a 
35% pileup fraction. We therefore extracted annular regions 
wherein we excised the « 20" radius central region. The ex- 
cised data accounted for approximately 1/3 of the detected 
events. The outer radii of our annular extraction regions were 
limited by the 1/4 sub-array used in our CCD readout mode, 
and thus were « 2'. We estimate that the extracted events had 
< 4% mean residual pileup fraction. All regions of the CCD 


are dominated by source counts, therefore we did not extract 
nor use any background spectra for the XIS observations. 

Events in the XIS detectors were read out from either 3x3 
or 2 x 2 pixel islands. The XIS 1 data were always in 3 x 3- 
mode, whereas the other detectors had mixtures of 3 x 3- 
and 2 x 2-mode. We created individual spectra and response 
files for each detector and data mode combination. Spec- 
tra were created for the time intervals described in §3, with 
the resulting exposure times also being listed in Table 1. Re- 
sponse matrices and effective area files were created with the 
xisrmfgen and xissimarfgen tools, respectively. 

Although for each observation the individual XIS spectra 
were fit separately, they were jointly grouped on a common 
grid such that they had a minimum combined signal-to-noise 
ratio of 8 in each energy bin (i.e., 64 total counts in each bin) 
and that the minimum number of channels per energy bin was 
at least the half width half maximum of the spectral resolu- 
tion 1 . To avoid regions of poorly understood response, we 
only considered spectral energy ranges 0.8-1. 72 keV, 1.88- 
2.19keV, and 2.37-7.5 keV 2 . 

For the PIN, we extracted spectra from the “cleaned” 
event files in the hxd/event__cl directories. The ap- 
propriate response and background files were downloaded 
from the High Energy Astrophysics Science Archive Re- 
search Center (HEASARC), specifically those from the 
pinxb_ver2 . 0_tuned directory. Good time intervals 
(GTI) were merged from our time selections (§3) and the 
combination of the GTI intervals from the source and back- 
ground event files. These intersected intervals were then 
used to extract the PIN source and background spectra. The 
source spectra exposure times were then corrected with the 
hxdtcorr tool. The PIN spectra were grouped to have a 
signal-to-noise ratio > 10 in each energy bin, and we consid- 
ered spectra between 12-70keV. 

The GSO spectra were created starting with the “unfil- 
tered” event files. These data were first reprocessed with 
the hxdtime, hxdpi, and hxdgrade tools, following 
the “ Suzaku ABC Guide”. The data were then filtered in 
xselect with the standard criteria from the HEASARC pro- 
vided gso_mkf . sel script. The background was down- 
loaded from the gsonxb_ver2 . 0 directory at HEASARC. 
GTI from the event file, the background file, and the time in- 
tervals were merged, and spectra were extracted from these 
times. Response files were then taken from the CALDB 
database, and exposure times were adjusted to agree with the 
spectra. The grouping of the GSO spectra is essentially fixed 
by the grouping of the background file, thus no rebinning was 
performed on these spectra. It was determined that the back- 
ground becomes prohibitive above 300 keV; therefore, we re- 
strict the GSO spectra to the 60-300 keV range. 

2.2. RXTE Analysis 

1 The half width half maximum (HWHM) was determined by using the 
spectral responses to create fake spectra of delta function lines at discrete 
enemies. These fake spectra were then fit without response matrices to de- 
termine HWHM. In practice, this meant that the original spectra with 4096 
channels each were grouped, using the I S I S group function, by a minimum 
of 6, 8, 12, 14, 16, 18, 20, 22 channels starting at 0.5, 1, 2, 3, 4, 5, 6, and 
7 keV, respectively. For the most part, the spectral binnings were dominated 
by the HWHM criteria rather than the signal-to-noise criterion. 

2 Outside of these energy ranges there are large disagreements among the 
fit residuals for the individual detectors and data modes. The excised regions 
between 1. 72-2.37 keV correspond to poorly calibrated Si and Ir features re- 
lated to the detectors and mirrors. 
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The RXTE data were prepared with tools from the HEA- 
SOFT v6.8 package and the most current calibration files as 
of 2010 January 10. We used standard filtering criteria for 
data from the Proportional Counter Array (PCA; Jahoda et al. 
1996). Specifically, we excluded data from within 30 minutes 
of South Atlantic Anomaly (SAA) passage, from whenever 
the target elevation above the limb of the earth was less than 
10°, and from whenever the electron ratio (a measure of the 
charged particle background) was greater than 0.15. We used 
the background models appropriate to bright data. 

We applied 0.5% systematic errors to all PCA channels, 
added in quadrature to the errors calculated from the data 
count rate. For all PCA fits we grouped the data starting at 
> 3 keV with the criteria that the signal-to-noise (after back- 
ground subtraction, but excluding systematic errors) in each 
bin had to be > 4.5. We restricted the noticed energy range to 
3-22 keV. 

We extracted data from the High Energy X-ray Timing Ex- 
periment (HEXTE; Rothschild et al. 1998) using the same 
criteria as for the PCA. HEXTE is comprised of two clus- 
ters, A and B. Prior to 2007, both clusters were in a rocking 
mode, with only one detector viewing the source at a given 
time, while the other detector conducted off-source back- 
ground measurements. For later dates (e.g., our three most 
recent observations), HEXTE A has been in a fixed on-source 
position, which required us to estimate cluster A backgrounds 
using the hextebackest tool. Furthermore, during our 
third observation there was clear evidence for a HEXTE A 
exposure anomaly 3 . We found evidence for such an anomaly 
(albeit less severe) in our second observation as well. We cor- 
rected the HEXTE A exposure time by iteratively fitting an 
exponentially cutoff powerlaw simultaneously to both cluster 
A and B data and varying the exposure of cluster A until a 
X 2 minimum was achieved using the same normalization for 
both detectors. 

We binned the HEXTE spectra to a common grid such that 
they had a minimum combined signal-to-noise of 8 in each en- 
ergy channel (even though both clusters were fit individually). 
We further restricted the noticed energy range to 1 8-200 keV. 
Additionally, when fitting the X-ray spectra we allowed the 
normalization of the HEXTE backgrounds to vary. (The best 
fit normalization constants were typically within < 10% of 
unity.) 

For both PCA and HEXTE spectra we further restricted the 
considered time intervals to those that were strictly simultane- 
ous with the Suzaku spectra. That is, the RXTE time intervals 
form a subset of the Suzaku time intervals. 

2.3. Chandra-HETG Analysis 

The High Energy Transmission Grating (HETG; Canizares 
et al. 2005) was inserted for our observation of Cyg X-l , with 
the data readout mode being Timed Exposure-Graded. The 
HETG is comprised of the High Energy Gratings (HEG), with 
coverage from a 0.7-8keV, and the Medium Energy Grat- 
ings (MEG), with coverage from « 0.4—8 keV. To minimize 
pileup in the gratings spectra, a 1/2 subarray was applied to 
the CCDs. Additionally, the observatory aimpoint was placed 
closer to the CCD readout. This configuration reduces the 
frame time to 1.741 sec, without any loss of the dispersed 
spectrum. 

We used CIAO v4 . 2 and CALDB v4.2.0 to extract 

3 http : //gsf c . nasa . gov/docs /xte/whatsnew/ 

newsarchive 2007. html/ 



Time (sec) 

Fig. 1. — Suzaku-XIS 1 lightcurves (16 sec integration bins) for observa- 
tion 4 in the 0.5-1. 5 keV band (filled diamonds, bottom curve) and the 3- 
9keV band (hollow diamonds, top curve, shifted upward by 50 cps). 

the data and create the spectral response files. The loca- 
tion of the center of the 0 th order image was determined us- 
ing the findzo.sl routine 4 , which provides » 0.1 pixel 
(^ 0.001 A, for MEG) accuracy. The data were reprocessed 
with pixel randomization turned off, but pha randomization 
left on. We applied the standard grade and bad pixel file fil- 
ters, but we did not destreak the data. 

Although the instrumental set up was designed to minimize 
pileup, it is still present in both the MEG and, to a lesser ex- 
tent, the HEG spectra. We incorporate the effects of pileup in 
our spectral fits using the simple_gpile2 model (Hanke 
et al. 2009) based upon the model originally described by 
Nowak et al. (2008). For the small amount of pileup present 
in these data, this model is extremely successful in describing 
its effect on the spectra. 

As for the PCA and HEXTE spectra, we restricted the con- 
sidered time intervals for the Chandra spectra to those that 
were strictly simultaneous with the Suzaku spectra. 

2.4. Estimated Radio fluxes 

Simultaneous Ryle radio fluxes are not available for the 
X-ray observations discussed in this work. Owing to the 
known stringent correlations between the X-ray and radio 
flux, specifically the correlations between the RXTE -ASM 
and Ryle fluxes (see Gleissner et al. 2004; Nowak et al. 2005), 
we use an estimated radio flux when applying the jet model 
to the spectra (§5.2). From our prior observations we average 
all radio data that occurred at times with ASM daily aver- 
ages that were within ±1 <j of the ASM daily average corre- 
sponding to these new observations. The average and stan- 
dard deviation of these radio data were used in the spectral 
fits, which yielded radio fluxes of 12.6 ± 3.3 mJy (2006-10- 
30), 10.8 ± 3.9 mJy (2007-04-30), 7.9±3.0mJy (2007-05- 
17), and 10.2±3.0mJy (2008-04-19). Each average was 
comprised of a minimum of ten measurements, and the ex- 
treme radio flux values within each sample never differed 
from the average by more than a factor of two. 

2.5. Data Plots 

4 http : //space .mit . edu/ASC/analysis/f indzo/ 
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Fig. 2. — Suzaku -XIS 1 color-color diagrams, made with 16 s bins. The colors are created from the count rate in three bands: 0.5-1.5keV, 1.5-3 keV, and 
3-9 keV. The lines are the best fit dust halo models, presuming that the source is obscured by both an interstellar neutral column as well as by a neutral column 
local to the source, with the latter ranging in value from 0-2 x 10 23 cm -2 . Filled symbols indicate the periods used for the spectral analyses described in this 
work. (See text.) 


Throughout this work we will present spectra that are “un- 
folded” from the detector response in a model-independent 
manner (see Nowak et al. 2005). Specifically, we define the 
unfolded photon flux, iw o id(/0, from a pulse height (PHA) 
bin, h 9 as 

p C(h)-B(h) 

Funfbld(A) f R(h,E)A(E)dE At ’ ( ) 

where C(h)-B(h ) are the background subtracted counts, 
R(h,E) and A(E ) are the detector response matrix and effec- 
tive area, and At is the exposure time. For multiple data 
sets coadded together for the plots (e.g., the Suzaku-XIS or 
the RXTE-HEXTE data), the numerator and denominator are 
each replaced with the sum of the same quantities from the 
detectors. Throughout we only sum data for the plots, not for 
the fits. 

Model spectra are “unfolded” in a similar manner. That is, 
predicted background subtracted model counts are first calcu- 
lated, and then unfolded as above. Thus, the model spectra 


exhibit the smearing induced by the detector response matrix. 
The one exception to this model component scheme is for the 
jet results presented in §5.2. There some of the model compo- 
nents are shown unsmeared by the detector response, as they 
are derived from internal calculations of the jet code. 

3. LIGHTCURVES 

By design, our fourth observation occurred near binary or- 
bital phase 0 (superior conjunction of the black hole), wherein 
we expected to be viewing Cygnus X-l through the wind of 
the secondary. By happenstance, all of our remaining Suzaku 
observations also occurred near superior conjunction, i.e., ob- 
servations 1-4 cover orbital phases 0.2-0.3, 0.8-0.0, 0.85- 
0.0, and 0.0-0.14, respectively. Although dipping events are 
known to be present during all orbital phases, during upper 
conjunction dipping is prevalent due to absorption by clumps 
in the wind (Bahicinska-Church et al. 2000). This behavior 
is clearly seen in the Suzaku lightcurves presented in Fig. 1. 
The soft X-ray lightcurve shows a deep, prolonged « 12 ksec 
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dip. A less pronounced, but still significant dip is seen in the 
hard X-ray lightcurve (Hanke et al. 2008). Numerous dips of 
varying levels are also seen on even shorter time scales. In 
order to extract the unabsorbed spectrum of the source, it is 
therefore necessary to screen the data for dipping and excise 
the dipping intervals from further analysis. 

In addition to dipping, spectral modeling also has to take 
into account the strong dust scattering halo in front of Cyg X- 
1 (Predehl & Schmitt 1995). As discussed by Xu, McCray 
& Kelley (1986, and references therein), a scattering halo is 
produced by the scattering of radiation from an X-ray source 
by a foreground dust cloud. Typical halo sizes are on the order 
of arcminutes, i.e., comparable to the size of the Suzaku PSF. 
The energy dependence of the scattering cross section leads 
to a halo angular size that is oc E~ l and a halo spectrum that 
is ex E~ 2 . Furthermore, as the scattered photons travel along 
a greater path length to reach the observer, they are delayed 
with respect to the directly observed photons by a factor that 
is proportional to the square of the angular radius from which 
the scattered photons are observed. The expected time delays 
for nearby galactic sources typically range from thousands to 
tens of thousands of seconds. 

Thus the Suzaku spectrum and lightcurve of Cyg X-l are 
comprised of two components: a directly observed compo- 
nent, subject to the local dipping events, and a time-delayed 
and time-averaged (due to integration over different angular 
radii) softer spectrum from the scattering halo. Most spec- 
tral analyses in the past have ignored the effects of scattering 
since to first order for optically thin scattering in a homoge- 
neous cloud, whatever radiation that is scattered out of the 
line of sight is scattered back in from larger (often spatially 
unresolved) radii. Thus, for spatially unresolved (on the size 
scale of the halo) and steady sources, dust can be ignored in 
simple analyses. 

If one assumes that the time-delayed and time-averaged 
spectrum is represented by the same model as the direct spec- 
trum, and if one further assumes that interstellar absorption 
is predominantly in the foreground of the dust halo, within 
ISIS one can write the spectral model for the spectra during 
the dips as: 

TBnew(l) * (1 - dust scat * (1 - TBnew(2)))* 

(cont inuum model) . (2) 

Here and throughout we shall use an updated version 5 * of 
the absorption model of Wilms, Allen & McCray (2000) to 
describe both the interstellar (model instance 1) and local 
(model instance 2) absorption. The dust scat model is a 
version of the XSPEC dust model that removes the opti- 
cally thin assumption of the latter model, but otherwise ig- 
nores multiple scatterings (F. Baganoff, priv. comm.). Note 
that when local absorption is absent (TBnew (2 ) — ► 1), the 
dust scattering term drops out of the model expression. 

Not all of the assumptions encompassed within the above 
model expression are necessarily realistic. Specifically, the 
assumptions that all interstellar absorption is foreground to 
the dust and that the delayed source spectrum is identical to 
direct spectrum are undoubtedly wrong in fine detail. How- 
ever, we find that the above expression provides a reasonable 
description of the color-color diagrams in Fig. 2. Specifically, 
we create a spectrum from bright phases of the lightcurve that 

5 http : //pulsar . sternwarte .uni-erlangen . de/wilms/ 
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TABLE 2 

Observational Fluxes in keV Energy Bands 


Date 

(yyyy-mm-dd) 

0.5-2 

2-10 

(10“ 8 

10-100 
erg cm -2 s -1 ) 

100-300 

LBoi a 

(LEdd) 

2006-10-30 

0.16 

0.78 

2.76 

1.59 

0.032 

2007-04-30 

0.07 

0.55 

2.07 

1.32 

0.026 

2007-05-17 

0.07 

0.38 

1.33 

0.87 

0.017 

2008-04-19 

0.14 

0.72 

2.62 

1.44 

0.035 


Note. — Fluxes are absorbed values, are normalized to the PC A spectral 
fit, and correspond to the brightest/least absorbed periods of the lightcurves. 
a Unabsorbed, isotropic luminosity in 0.01-800 keV band, expressed as a 
fraction of Eddington luminosity, as determined from the thermal Comp- 
tonization fits of §5.1 and assuming a distance of 2.3 kpc and a black hole 
mass of 10Mq. 

inhabit a locus of points in the upper right hand comer of the 
color-color diagrams (i.e., the filled symbols of Fig. 2). We 
fit a simple spectral model to these selected data: an absorbed 
disk plus powerlaw plus broad and narrow gaussian lines near 
6.4 keV. We then apply this model in eq. (2), and create color- 
color curves by varying the TBnew (2 ) component from 0- 
2 x 10 23 cm -2 . These curves are then fit to the color-color 
diagrams using the optical depth of the dust scattering halo as 
the single free parameter. Our best fits yield r = 0.24—0.34, 
and the results are shown in Fig. 2. 

Overall, these curves describe the behavior of the color- 
color diagrams reasonably well. The dust halo represents 
» 20%-30% (i.e., on the order of the scattering optical depth) 
of the soft X-ray (^ 2 keV) flux being “uncovered” during 
the dips. Given the much narrower Chandra PSF, one would 
expect comparable Chandra color-color diagrams to show a 
far smaller uncovered fraction if the dust halo interpretation 
is correct. This is indeed found to be true, with Chandra 
showing only an « 2% uncovered fraction (Hanke et al. 2008, 
Hanke et al., in prep.). However, it is possible that the very 
low uncovered fraction required by the Chandra data might, 
in fact, be partly intrinsic to the local dipping at the source 
(Hanke et al. in prep.). Note that the theoretical curves show 
a greater degree of curvature than the observed, nearly linear, 
evolution toward the lower left hand comer of the diagram. 
We hypothesize that this discrepancy is due to our assumption 
in eq. (2) of purely neutral absorption in the dips. A concur- 
rent increase in the optical depth of ionized absorption could 
serve to alter the theoretical curves in the appropriate manner. 
We return to this concept of an ionized absorber component 
in §4.2. 

4. SIMPLE SPECTRAL MODELS 

The lightcurves and color-color diagrams described above 
were used to create spectra for the time periods corresponding 
to the bright phases of the lightcurve represented in the upper 
right comer of the color-color diagrams (Fig. 2). The selec- 
tions that we made are highlighted as filled symbols within 
these figures. The resulting exposure times are given in Ta- 
ble 1 . With the exception of the RXTE - only broken powerlaw 
spectral fits described immediately below, throughout the rest 
of this work we shall only describe spectra from these bright 
phase, minimally locally absorbed time periods. 

4.1. Broken and Cutoff Powerlaw Descriptions 

As discussed by Wilms et al. (2006), nearly all RXTE- 
PCA and HEXTE spectra of Cyg X-l can be described by 
an extremely simple phenomenological model: an absorbed, 
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TABLE 3 

CONSTANT* HIGHECUT* POWERLAW PARAMETERS FOR HARD X-RAY SPECTRA FITS 


Date 

(yyyy-mm-dd) 

r 

-4pl 

£eut 

(keV) 

£fold 

(keV) 

Qjso 

Chexte a 

Chexte b 

onorm 
° HEXTE A 

Dnorm 
U HEXTE B 

X 2 /D oF 

2006- 10-30 

2007- 04-30 

2007- 05-17 

2008- 04-19 

i 41 + 0.02 

j 42+8:81 

l 47+8:81 

A>t/ -Q.Q 2 
1 aQ+0.01 

A - jy -o.oi 

1 07 +005 

q 73+8:81 
u * /J -Q.Q3 

u - ou -Q.Q3 

Q 03 + 0.02 

u - yj -0.03 

19-618 

24.6 j ^ 
23-5 If 
20.4+] j 

183 + ? 

218+' 

252 + ?l 

164*! 

1 007* 0008 
095« 8 

q 95+8:81 

Q 01 +0.01 
uyA -0.01 

1 -O^of 

n o-i 7+0.003 
U.5J /_o 003 

0 OO7+0-005 
-0-004 
1 07+0.02 
AU/ -0.01 

0 842 + 0 004 
U.5^Z_0 004 

1 1 7 +O.OI 

1 03+8:88 
A *^ -0.06 

1 065+ 0008 

1 02 40 02 

0 99+8:84 
Uyy -0.06 

0 QQ+ 0.01 

uyy -o.oi 

785.4/449 

966.2/315 

301.6/123 

648.8/447 


Note. — Spectra are normalized to HXD-PIN data, ^cut, £foid are the cutoff and folding energy of the highecut model. A p i is the power law normalization 
in photons keV -1 cm -2 s -1 at 1 keV. C are the fit constants for the detectors other than the HXD-PIN. The normalizations of the HEXTE backgrounds were also 
allowed to be free parameters, and are given by the fi nonn values. 


exponentially cutoff, broken powerlaw plus a broad, gaus- 
sian line. The photon indices of the broken powerlaw, which 
we label as T\ and T 2 for the soft 3-10keV) and hard 
(<; lOkeV) X-rays, show a correlation in the sense that the 
amplitude of the break between the two, T \ -T 2 , increases 
with higher (i.e., softer) values of r 1# This T-AT correla- 
tion is a very phenomenological description of the spectral 
correlation that elsewhere has been described as a ‘T-fi/27r” 
(hardness/reflection traction) correlation (Zdziarski, Lubinski 
& Smith 1999). Following Remillard & McClintock (2006), 
Wilms et al. (2006) classified spectra with T i < 2.2 as “hard 
state” spectra. 

In order to compare the spectra discussed here to the work 
of Wilms et al. (2006), we take the 9 RXTE spectra 6 * and fit 
them with exponentially cutoff broken powerlaws. As in the 
work of Wilms et al. (2006), these simple phenomenological 
models are excellent descriptions of these data. This is even 
true for the spectra extracted from times of moderate and deep 
dips. These latter spectra have slightly reduced values of Tj 
and increased values of the fitted neutral column, Nr (Wilms 
et al. 2006 show that orbital phase-dependence of the column 
is discernible in our RXTE observations of Cyg X-l), but oth- 
erwise appear similar to non-dip phase spectra. This is not 
entirely surprising given that RXTE is primarily sensitive at 
energies <; 3 keV, and therefore does not measure the spectral 
regime most strongly affected by the dips. 

We show the resulting photon indices of these spectra 
along side the results from Wilms et al. (2006) in Fig. 3. 
Our prior study covered seven years and 202 RXTE spec- 
tra. These newer observations overlap the historically hard- 
est states measured in that campaign. Specifically, for these 
new observations T i ranged from 1.63-1.71, which is to be 
compared to the T\ = 1.65 minimum found in Wilms et al. 
(2006). In some cases, these observations also correspond to 
historically faint hard states of Cyg X-l. Defining the “hard 
state” as Ti < 2.2, the lowest 2-100keV flux in the hard state 
reported by Wilms et al. (2006) was 2.46 x 10“ 8 ergcm _2 s _1 , 
while the highest 2-100keV hard state flux reported in that 
work was 4.50 x 10 _8 ergcm _2 s -1 . As shown in Table 2, the 
2-100keV fluxes spanned by the bright phase, least absorbed 
observations range from 1 .7 1-3 .54x1 0 -8 erg cm -2 s -1 . 

That these spectra are among the faintest and hardest ever 
observed for Cyg X-l is important to note. In some models 
of the black hole hard state, the hardest spectra are expected 
to have the smallest observed reflection fraction and the disk 

6 These 9 spectra consist of three bright phase spectra, and six spectra 

taken from deeper dipping periods. Observation 1 shows no strong dips, and 

observation 3 does not have RXTE data strictly simultaneous with the bright 
phases of the Suzaku spectra. 



Fig. 3. — Correlation of the soft X-ray (Ti) and hard X-ray (T 2 ) photon 
indices from exponentially cutoff, broken powerlaw fits to the RXTE Cyg X- 
1 data presented by Wilms et al. (2006). Triangles represent the hard states 
(Ti < 2.2), and circles represent soft states (T 2 > 2.2). The solid diamonds 
represent exponentially cutoff, broken powerlaw fits to the RXTE data dis- 
cussed in this work. There are spectra from four observations that have been 
further subdivided into up into times covering different portions of the color- 
color diagrams of Fig. 2, for a total of 9 RXTE spectra. 

inner radius should be at a maximum (Zdziarski, Lubinski & 
Smith 1999). However, it is also important to note that other 
black hole systems can show substantially fainter spectra rela- 
tive to Eddington luminosity. Furthermore, as we see here and 
as has been noted elsewhere, the range of hard state fluxes ob- 
served in Cyg X-l spans only a factor of two, and it is unlikely 
that the bolometric flux spans a much greater factor than this 
even when including observations from the soft state (Wilms 
et al. 2006). 

As a further comparison to the previously observed hard 
states, we consider the value of the fitted folding energy. For 
the hard states discussed by Wilms et al. (2006), the fold- 
ing energy ranges from « 125-255 keV. Other BHC, e.g., 
GX 339-4, have shown a wider range of hard state fold- 
ing energies, i.e., ^ 50-300 keV (Wilms et al. 1999; Nowak, 
Wilms & Dove 2002; Nowak et al. 2005; Motta, Belloni & 
Homan 2009). In our prior studies we did not consider spec- 
tra above 125keV. Here, however, with the inclusion of the 
Suzaku -GSO spectra we can now consider energy ranges up 
to 300 keV. For all four of our observations, an exponential 
rollover is clearly detected. An example is given in Fig. 4, 
which also gives a general indication of the consistency of the 
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Fig. 4. — Suzaku-HXD PIN (brown squares) and GSO (magenta trian- 
gles) spectra, and RXTE-UEXTE (green diamonds) spectra from the least 
absorbed periods of observation 4 fit with an exponentially cutoff powerlaw. 

cross-correlation among the hard X-ray detectors. The values 
of the folding energies, along with cross-normalization con- 
stants, are given in Table 3. The folding energies range from 
164— 252 keV, nearly the full span encompassed by our prior 
hard state observations. This is not necessarily surprising. 
Whereas Wilms et al. (2006) had observed a general trend for 
the cutoff to increase to higher energies with harder spectra, 
there was a large degree of scatter about this trend, consistent 
with these four observations. 

4.2. Ionized Absorption Models 

Previous Chandra -HETG observations of Cyg X-l have re- 
vealed evidence of the ionized, focused wind from the donor 
star (Miller et al. 2005; Hanke et al. 2009). The effects of 
this wind upon the spectrum must be accounted for in any de- 
tailed analysis of Cyg X-l. While ionized line absorption is 
very pronounced during the deepest part of the absorption dips 
(Fig. 1), we stress that such absorption is also significant dur- 
ing the bright phase (non-dip) spectra discussed in this work. 

We begin by considering the joint 0.8-7.5keV Suzaku 
and Chandra -HETG spectra from the 2008 April 19 ob- 
servations. As the MEG spectra are more heavily af- 
fected by photon pile-up, we fitted only the HEG spec- 
tra. We used the same absorbed disk plus powerlaw 
(TBnew* (diskbb+powerlaw) ) continuum model to fit 
both the Suzaku and HETG spectra. We also added to 
both spectral models a broad Fe Ko; line (modeled with 
diskline; Fabian et al. 1989) and a narrow Fe Ka line 
modeled with a gaussian, and a 6.7 keV (Fe xxv) absorption 
line. Following Hanke et al. (2009), we describe the line ab- 
sorption as a series of gaussian absorption lines representing 
a variety of ionized species — predominantly H- and He-like 
lines from elements ranging from O to Fe, and L-shell transi- 
tions of Fe. This initial line list contained nearly 100 line tran- 
sitions parameterized by wavelength (which we constrained to 
fall within 1500 km s" 1 of their rest wavelengths), full width 
half maximum (constrained to be between 5^10 mA), and line 
equivalent width. Based upon fits to the joint Suzaku-HETG 
spectra, this initial list was reduced to the 55 significant lines 
present in the spectra. 

The above joint model received two modifications for solely 



Energy (keV) 

Fig. 5. — Suzaku (blue diamonds) and Chandra -HETG (grey histogram; 
d= 1 st order HEG only) data from observation 4. The model consists of ab- 
sorbed disk and powerlaw components, broad and narrow gaussian lines in 
the Fe Ka region, and absorption lines from a variety of highly ionized 
species, including from Fe xxv at « 6.7 keV. Additionally, we modify the 
HETG spectra by dust scattering and a pileup term (see text). The Suzaku 
and HETG spectra have not been renormalized with respect to one another. 
The residuals shown here are for this joint model fit with all broad and narrow 
line normalizations (both emission and absorption) set to zero. 

the HETG spectra. The simple_gpile2 model (see Hanke 
et al. 2009) was applied to describe the effects of pile up in 
the gratings spectra. Owing to this spectrum being fainter 
and harder than the Cyg X-l spectrum described in Hanke 
et al. (2009), these spectra required an even smaller pileup 
correction. Additionally, we applied the dust scattering model 
(dustscat) to the HETG spectrum. As described in §3, 
dust scattering represents a loss term for the high spatial res- 
olution Chandra spectra. This component, albeit with a time 
delay of thousands to tens of thousands seconds, scatters back 
into the Suzaku spectrum from arcsecond to arcminute angu- 
lar scales. To the extent that one can ignore the time evolution 
of the spectrum, the dust scattering term can be ignored in the 
Suzaku spectra. Lacking any detailed information of the aver- 
age spectrum prior to the start of our observations, we do not 
apply any dust scattering correction to the Suzaku spectra. 

To be explicit, we apply to the Suzaku data a model of the 
form: 

TBnew* lines* 

(di skbb + powerlaw 4- di ski ine + gaus s i an) (3) 

while we apply to the HETG data a model of the form: 

simple jrpile2(g) (TBnew* lines * dustscat* 

(di skbb + powerlaw + di ski ine + gaus si an)) , (4) 

where lines represents the ionized line absorption. 
Throughout the rest of this work, we shall ignore as we are 
doing here the effects of dust scattering on the RXTE and 
Suzaku data. Again, this implicitly assumes that the bright 
phase spectrum has been steady over the time span of thou- 
sands to tens of thousands of seconds. This assumption is 
unlikely to be valid in detail; however, we lack the data to 
employ any more sophisticated assumptions. 

This simple model describes the spectra well with a dust 
halo optical depth of 0.25 =L 0.03 (90% confidence level), i.e., 
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Fig. 6. — The HEG data from a joint fit of the Suzaku and Chandra data 
with a simple continuum model (disk, powerlaw, broad and narrow Fe re- 
gion lines, neutral ISM absorption) plus two ionized absorber components 
described with the warmabs model. The residuals include all of these mod- 
els components; the remaining line residuals are features not described by the 
warmabs components. 

consistent with the fits to the color-color diagrams presented 
in §3. This spectral fit is shown in Fig. 5. The important 
point to note here is that the ionized absorption is extremely 
statistically significant in the Suzaku spectra, even though the 
individual lines are not resolved. Inclusion of these lines is vi- 
tal for obtaining a description of the soft end of the spectrum, 
and thus we include ionized absorption in all of our model fits. 
Being that we do not have HETG spectra simultaneous with 
our other three Suzaku spectra, we use the ionized absorption 
fit presented in Fig. 5. Specifically, for all subsequent spec- 
tral fits described in this work, we freeze the line positions, 
widths, and relative equivalent widths to the values found for 
this joint fit. The line equivalent widths are tied together via 
a single normalization constant, and this normalization con- 
stant becomes the sole fit parameter describing the ionized 
absorption. The model, however, still consists of 55 individ- 
ual absorption lines. As described below, this approach works 
well in all of our spectral fits. 

We explored whether a more sophisticated ionization model 
could be applied to our data. Rather than fit individual lines 
to the joint Suzaku and Chandra data, we instead used the 
warmabs model (Kallman & Bautista 2001) to describe the 
ionized absorption. A model that includes two ionized com- 
ponents — one with a column of (7 ± 4) x 10 21 cm -2 and ion- 
ization parameter log 10 £ = 2.77 =L 0.05 and one with a col- 
umn of (3 =t 2) x 10 21 cm -2 and ionization parameter log 10 £ = 
1.93 ±0.07 (90% confidence levels), and both with densi- 
ties of 10 10 cm“ 3 — describes the data reasonably well. As 
shown in Fig. 6, however, a number of prominent line fea- 
tures are not well-fit in the Chandra data. Furthermore, aside 
from being an extremely computationally expensive model to 
run, when fit to the Suzaku data alone the warmabs mod- 
els would tend to gravitate toward low ionization parameter 
values that, although capable of mimicking some of the con- 
tinuum features of the Suzaku data, did not include the high 
ionization line features that were clearly present in the Chan- 
dra data. For these reasons, we choose the empirical approach 
outlined above. 



Energy (keV) 

Fig. 7. — The summed unfolded Suzaku spectra from unabsorbed periods 
from observations 2-4, shown in the Fe Ka line region. These data have been 
fit with a model consisting of an absorbed powerlaw, narrow Fe Ka emis- 
sion line, narrow Fe xxv absorption line, and a relativistically broadened 
diskline model. (See text for fitting procedures.) The individual model 
components (powerlaw, with added diskline, and separately with added nar- 
row emission and absorption lines) are shown overlain on the spectra, while 
the residuals are shown with and without the emission and absorption line 
components. 

Several things need to be borne in mind when considering 
our results: whereas we allow for an overall normalization 
change in our ionized absorber model, we do not allow for 
changes of the ionization state. As Fig. 5 shows, with the 
exception of a weak Fe xxv line, there are essentially no 
ionized lines with energies ^ 3 keV. Thus it is unlikely that 
this simple empirical approach will directly impact, for exam- 
ple, broad Fe line studies. However, in as much as ionization 
changes could affect our soft X-ray continuum fits, and these 
continuum fits in turn affect our estimate of the breadth and 
strength of the red wing of any broad line tail (see §5.3), po- 
tential ionization changes in this absorber must be considered 
an additional source of systematic uncertainty in the results 
discussed below. 

4.3. Composite Relativistic Line 

The fit described above consisted of both broad and narrow 
features in the Fe Ka region. Here we examine this spec- 
tral region in a somewhat phenomenological manner to gauge 
the interplay between these narrow and broad features. We 
wish to consider the average line profile from roughly simi- 
lar Cyg X-l spectra. The three most recent Suzaku observa- 
tions occurred closest to orbital phase 0 and show very similar 
color-color diagrams. As discussed below, these three spectra 
are fit with comparable ionized absorption normalization con- 
stants. We therefore consider these spectra jointly and ignore 
the first Suzaku observation. We fit these spectra with a model 
consisting of a powerlaw, narrow emission near 6.4 keV and 
narrow absorption near 6.7 keV, plus a relativistically broad- 
ened Fe line (again using the diskline model). The rela- 
tivistically broadened line is characterized by a broad tail ex- 
tending redward from the line rest frame energy, and a sharper 
peak (due to Doppler boosting) blueward of the line rest frame 
energy. 

For purposes of this discussion, we allow the powerlaw to 
assume independent parameter values for each of the Suzaku 
observations. The energies and widths of the narrow lines are 
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Fig. 8. — Suzaku (blue line), Chandra-HEG (grey line), and RXTE - PCA 
(purple diamonds) residuals from fitting an absorbed powerlaw to the 3.5- 
4.5 keV and 7-7.5 keV regions simultaneously. The Chandra-HEG data have 
been binned to match approximately the binning of the Suzaku data. 

tied together for all three observations, but the line strengths 
are fitted individually. The broad line has an energy fixed 
to 6.4 keV and an emissivity index fixed to /3 = -3 (i.e., line 
emissivity is oc R~ 3 9 where R is the emission radius within 
the disk), and the disk inclination is fixed at 35° (the value 
adopted by Herrero et al. 1995, from the middle of the range 
of suggested inclinations). The inner radius of the broad line 
emission is tied together for all three observations, but the line 
strengths are left independent. The x 2 of the fit is determined 
from the sum of \ 2 values from the individual observations; 
however, in Fig. 7 we show the spectra, fitted models, and 
residuals from the combined data. We fit the spectra only in 
the 4.5-7. 5 keV bandpass. 

These joint fits strongly require the presence of a broad 
line. The broad line is significantly detected in each individ- 
ual XIS detector from each individual observation. Under the 
assumptions of the continuum model applied here, the broad 
line equivalent width is approximately three times that of the 
narrow line equivalent width. The best fit inner radius for 
the broad line emission is 10. 6^9 GM/c 2 (90% confidence 
level), i.e., close to the innermost stable circular orbit for a 
Schwarzschild black hole. As we will discuss further in §5.1, 
the specific value for such an inner radius is dependent upon 
the fitted continuum. Nevertheless, a broad excess is required 
redward of the Fe Ka line. 

In describing the residuals in this region with a broad line, 
we see that the Fe xxv absorption line occurs near the blue 
wing peak of the relativistically broadened line (Fig. 7). Both 
the narrow emission and absorption lines in the Fe region 
are required in the HETG spectra when considered by them- 
selves (see also Hanke et al. 2009), and thus they cannot be 
ignored in the Suzaku spectra. The HETG spectra, however, 
are not particularly well-suited for describing any broad com- 
ponent of the Fe line. (See Fig. 8 and the discussion of Hanke 
et al. 2009 concerning broad line fits to our prior joint RXTE- 
Chandra-HETG observations of Cyg X-l.) 

On the other hand, if we compare the RXTE spectra to both 
the Suzaku and HETG spectra, we find that the RXTE spec- 
tra alone require an Fe line equivalent width that cannot be 
accommodated with solely the narrow line component from 



Fig. 9. — Possible X-ray emission geometry of the Cyg X-l system. Part 
of the soft X-rays likely come from an accretion disk, which may he trun- 
cated at its inner edge hy a hot corona. Two possible (of many suggested) 
geometries are shown for this corona: a torus, and a quasi-spherical cloud. In 
some models, the disk can reform in the inner regions of the accretion flow 
(Mayer & Pringle 2007). The corona will Compton upscatter a fraction of the 
disk component(s) soft X-rays to hard X-rays. Whether or not the corona en- 
velopes the disk(s) determines whether or not an additional, unComptonized 
disk component should he present in the spectra. (If the disk reforms on the 
inner edge of a toroidal corona, for example, a fraction of its soft X-rays 
will not intersect the corona.) Hard X-rays will be reprocessed into reflec- 
tion/fluorescent line components, predominantly hy the outer disk. Finally, 
the jet might contribute to the observed X-ray emission via synchrotron and 
synchrotron self-Compton emission. 

the Suzaku-HETG fit. We highlight this fact in Fig. 8 , where 
we show the residuals from a powerlaw fit in the 3. 5-4.5 keV 
and 7-7.5 keV region to the joint RXTE -Suzaku -HETG spec- 
tra. In this figure, we have binned the HETG spectra to match 
approximately the binning of the Suzaku spectra. There is 
overall good agreement between the RXTE and Suzaku resid- 
uals; both require a broad red wing in the Fe line region. The 
narrow Fe Ka emission and Fe xxv absorption are clearly 
seen in the HETG spectra; however, when viewed by itself, 
any broad Fe line component in the HETG spectra could be 
subsumed via a slight powerlaw slope change. We ascribe 
this latter fact to remaining calibration issues in the HETG 
detector. Based upon these results, we include broad and nar- 
row Fe Ka emission and narrow Fe xxv absorption in all the 
subsequent fits presented in this work. As we shall discuss 
below, however, the implied broad line parameters are depen- 
dent upon the assumed continuum model. 

5. BROAD BAND MODELS 

We now turn to a discussion of the 0.8-300 keV, joint 
Suzaku-RXTE spectra for the bright phases of our four 
Cyg X-l observations. Here we consider three different 
Comptonization models, two of which are discussed in detail, 
and a jet model. (“Toy geometries” for some of these situa- 
tions are presented in Fig. 9.) The latter model also contains 
Comptonization components that represent both Comptoniza- 
tion of disk photons as well as synchrotron self-Compton 
(SSC). The non-jet models discussed below all rely upon the 
eqpair model (Coppi 1999, 2005) to describe the fitted 
Comptonization components. 

5. 1 . Comptonization Descriptions 

The eqpair model allows for Comptonization by a coro- 
nal electron population that has both thermal and non-thermal 
energy distributions. The latter distribution is governed by a 
parameterized powerlaw, which for purposes of the fits de- 
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Energy (keV) Energy (keV) 

Fig. 10. — Unfolded Suzaku and RXTE spectra, plus model components and fit residuals from the high kT see( j, thermal corona eqpair fits. Suzaku-XIS 
spectra from individual detectors/data modes and RXTE-HEXTE spectra from each cluster were summed for the figures, but not for the fits. Here, and throughout 
the remainder of this paper, XIS data are filled blue diamonds, PC A data are hollow purple circles, PIN data are hollow green diamonds, HEXTE data are brown 
squares, and GSO data are hollow magenta triangles. The following model components are shown: the summed model with the reflection component excluded, 
the reflection component by itself, the broad and narrow Fe Kq line components, the seed photons without Comptonization applied, the additional disk component 
with peak temperature tied to the seed photons, and the additional low temperature disk component. Neutral and ionized absorption are applied to all shown 
model components. 


scribed here we left at the default eqpair values. Specifi- 
cally, the electron phase space density follows a distribution 
oc E~ 2 , where E is the electron energy, between energies of 
1 .3 eV-1 MeV. The eqpair model rather than being parameter- 
ized by a coronal electron scattering optical depth and tem- 
perature is instead parameterized by a seed electron optical 
depth (pair production within the corona, accounted for in the 
code, can yield higher net scattering optical depths) and coro- 
nal compactness parameters, t (proportional to component lu- 
minosity or power divided by radius). The latter parameters 
are divided into two relative compactness parameters: 4/4, 
the coronal compactness divided by the seed photon compact- 
ness, and 4t/4, the compactness of the non-thermal electron 
population divided by the total electron population compact- 
ness. Over a wide range of parameters, the hardness of the 
spectrum and the temperature of the corona increases with in- 
creasing compactness (although pair production can modify 


that effect). 

The seed photon spectrum can be set to either a black- 
body or a disk spectrum (essentially the diskpn model 
from XSPEC; Gierlinski et al. 1999). We choose the latter 
seed photon spectrum, with the peak temperature of this disk 
(£7di s ki) becoming a fit parameter. We also allow for an un- 
scattered disk component by including an additional di s kbb 
component with peak temperature fixed to the peak temper- 
ature of the seed photons. Our choice for the normalization 
of this component depends upon the specific Comptonization 
model under consideration (see below). 

To account for soft excesses (i.e., below ^1.5 keV) in the 
spectra, we sometimes include an additional unscattered disk 
component with freely varying normalization and temperature 
(^disk 2 , *7 disk2 ). We searched for, and easily found, solutions 
where the peak temperature of this unscattered disk com- 
ponent was below « 300 eV. The eqpair model includes 
reflection from an ionized disk (i.e., a modification of the 
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Fig. 11. — Same as Fig. 10, but for the non-thermal corona, low kT^ eqpair model fits. The following model components are shown: the summed 
model with the reflection component excluded, the reflection component by itself, the broad and narrow Fe Ka line components, the seed photons without 
Comptonization applied, and the additional disk component with peak temperature tied to the seed photon peak temperature. Neutral and ionized absorption are 
applied to all shown model components. 


pexriv model; Done et al. 1992) which is smeared by rel- 
ativistic distortions with the emissivity profile of a Shakura- 
Sunyaev type disk (Shakura & Sunyaev 1973). A fittable pa- 
rameter is the inner radius of this disk smearing profile (see 
below). We also fix the disk (i.e., reflector) temperature to 
10 6 K, the reflector inclination to 35°, and allow the reflection 
fraction and reflector ionization parameter (£) to be fit 

parameters. The latter is limited to values < 1000. 

Line absorption is described by the parameterized (i.e., sin- 
gle normalization constant) model discussed in §4.2. A nar- 
row gaussian is added at 6.399 keV with its width frozen 
to 0.01 keV, but with a freely variable normalization. A 
diskline component (essentially the profile for a disk 
around a Schwarzschild black hole; Fabian et al. 1989) is 
added to describe the broad line. Its energy is fixed to 6.4 keV, 
the disk line emissivity index is set to j3 = -3, and the disk in- 
clination is set to 35°. Typically one allows the inner emission 
radius to be variable (see §4.3). Here we tie this inner radius 
to that of the reflector in the eqpair model. The line nor- 
malization, however, is allowed to freely vary. 


A number of authors have considered “sphere+disk” Comp- 
tonization models wherein an inner quasi-spherical corona is 
encircled by a cool, geometrically thin disk that typically has 
a peak temperature of « 200 eV (Gierlinski et al. 1997; Dove 
et al. 1997). We were unable to find any solutions for coronae 
with thermal electron distributions (i.e., £ n t/£h frozen at 10“ 3 ) 
that allowed for such a low seed photon temperature. The only 
such thermal corona solutions that we found required seed 
photon temperatures of 0.8-1 keV. Fit parameters for these so- 
lutions are presented in Table 4, and spectral fits are shown in 
Fig. 10. 

Concepts for the potential geometry represented by this 
spectral fit are shown in Fig. 9. The high temperature seed 
photons could be a recondensed inner disk as envisioned by 
Mayer & Pringle (2007), while the low temperature soft ex- 
cess, the broad line, and the reflection component could em- 
anate from an outer disk. (The narrow line likely arises from 
fluorescence from the secondary and/or spatially extended gas 
surrounding the system; Torrejon et al. 2010.) A more physi- 
cally self-consistent model would allow for an additional seed 
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Energy (keV) 

Fig. 12. — Same as Fig. 10, but only showing observation 4 for the mul- 
tiple eqpair component fit (see text). The following model components 
are shown: the two individual eqpair components with their reflection com- 
ponents excluded, the reflection component of the dominant eqpair spec- 
trum, the broad and narrow Fe Ka line components, and the unscattered disk 
component. 

photon contribution from this outer disk, and possibly a broad 
line component from the inner, hot seed photon emission re- 
gion. However, it is interesting to note that these solutions re- 
quired very little additional unscattered, high temperature disk 
component. (Table 4 contains low- sometimes zero- values of 
^4 disk i . Both the initial seed photon distribution and the associ- 
ated high temperature, unscattered disk spectrum are shown in 
Fig. 10; the former is much more significant.) The amplitude 
of the soft seed photon contribution and the weakness of an 
additional unscattered high temperature disk would indicate 
a ‘recondensed’ disk region of only modest extent, predom- 
inantly within the confines of a more extended corona, from 
which we would expect little iron line contribution. 

On the other hand, these spectral fits do require an addi- 
tional, unscattered low temperature disk component. The val- 
ues ofr4disk2 presented in Table 4 correspond to inner disk radii 
of « 2-10 GM/c 2 , given a distance of 2.3 kpc, a black hole 
mass of 1OM 0 , and an inclination of 35°. This is roughly 
consistent with the inner radius of the fitted broad Fe line 
and relativistically smeared reflector, for which we find val- 
ues ranging from R\ n = 6-1 8 GM /c 2 . (The inner emission ra- 
dius of the line appears to be uncorrelated with the inner ra- 
dius implied by the normalization of the low temperature disk 
component.) These radii values are roughly consistent with 
values close to the marginally stable orbit of a Schwarzschild 
black hole. 

The broad line amplitude is between 6-22 times the ampli- 
tude of the narrow line component. Fitted reflection fractions 
are Q/ln & 0.2. 

Comparable high seed photon temperature fits were found 
by Wilms et al. (2006) when considering solely RXTE spec- 
tra. The major difference between those fits and the ones dis- 
cussed here is that our prior fits also required a significant con- 
tribution from an unscattered, high temperature disk. It is un- 
clear whether that difference is due to the inclusion of Suzaku 
data that extend to lower energies, the inclusion here of ion- 
ized absorption, or the fact that these spectra are for the most 
part fainter and harder than those discussed by Wilms et al. 
(2006). The coronal compactnesses found here are compara- 


ble to the maximum values from our previous studies, while 
the seed optical depths are slightly lower (by « 0. 3^0.5 for 
comparable observations). These values of seed optical depth 
and compactness correspond to coronal temperatures of 85- 
1 15 keV and total optical depths of 1.1 -1.4. The fitted reflec- 
tion fractions are slightly higher here (by « 0.05) than those 
fits with comparable 4/4 from Wilms et al. (2006). 

We have found a set of Comptonization model solutions 
that do allow for a low seed photon temperature. For these 
fits, we include only one extra unscattered disk component 
(with its peak temperature tied to the peak seed photon tem- 
perature), and we further tie the inner disk radius to the inner 
radii of the eqpair reflector and the diskline emissivity. 
That is, we are explicitly mimicking the “sphere+disk” geom- 
etry (albeit not in a completely self-consistent manner) shown 
in Fig. 9. If we then relax the assumption of a purely thermal 
coronal electron distribution, good descriptions of the spec- 
tra are found with low seed photon temperatures. The ratio 
of non-thermal to total coronal compactness then falls in the 
range 4t/4 = 0.01-0.82. That is, up to 82% of the electron 
energy (presuming that the thermal and non-thermal electrons 
are cospatial) resides in a powerlaw distribution that extends 
to 1 MeV. Parameters for this model are presented in Table 5 
and spectra are shown in Fig. 1 1 . 

The quality of these fits is similar to or perhaps slightly 
better than that for the thermal corona fits. The fitted com- 
pactness parameters for the non-thermal coronae are approxi- 
mately double those for the thermal fits. (Similar results were 
found by Ibragimov et al. 2005 when comparing thermal and 
non-thermal eqpair fits.) At first glance, the optical depths 
seem smaller; however, owing to the large coronal compact- 
nesses pair production is significant and the net optical depths 
range from 0.7-1 .6, while the thermal electrons have temper- 
atures that range from 55-160keV. Reflection fractions are 
slightly larger (by up to 0.09) for these fits. The broad line in- 
ner radius in some cases has increased, and now ranges from 
» 6-32 GM/c 2 . Additionally, the amplitude of the broad line 
is now only 3-13 times that of the narrow line. This differ- 
ence compared to the thermal corona model is due to two fac- 
tors: the narrow line amplitude is slightly increased, while the 
broad line amplitude in two cases is nearly halved. 

We briefly consider one other thermal Comptonization so- 
lution that allows for a low seed photon temperature. (No “toy 
geometry” for this model is shown in Fig. 9.) It is a model 
similar to that considered by Frontera et al. (2001) when mod- 
eling BcppoSAX spectra of Cyg X-l and to that considered 
by Makishima et al. (2008) when modeling Suzaku spectra of 
Cyg X-l. In these models two eqpair components are in- 
cluded. They share the same seed photon temperature (tied to 
the peak temperature of an additional unscattered disk com- 
ponent, as before), and they share the same reflection parame- 
ters; however, they have independent compactness and optical 
depths. An example of such a model fit is shown in Fig. 12. 
These fits are generally successful, but no more so than the 
previous two models. Furthermore, owing to the larger degree 
of freedom given the two semi-independent coronal compo- 
nents, they are fraught with local minima representing qualita- 
tively different relative contributions of the two Comptoniza- 
tion components. We do not consider these more complex 
models farther in this work, other than to point out that such 
double corona solutions as discussed by Makishima et al. 
(2008) and Frontera et al. (2001) can allow for both low seed 
photon temperature and purely thermal coronae. 
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Fig. 13. — The jet model fit to the estimated radio flux and measured X-ray 
spectra from observation 4. Individual model components are shown with 
both neutral and ionized absorption applied. We show the summed model, 
as well as the the jet synchrotron component, the jet SSC component, the 
Compton corona component, and the disk component. 

5.2. Jet Descriptions 

The Compton corona models described above provide a 
good description of the 0.8-300 keV spectra. These models, 
however, are not unique. Even within the restricted class of 
coronal models, we have shown three qualitatively different 
solutions that yield comparable fits. Furthermore, none of 
the Compton corona models give a self-consistent description 
of the correlated radio spectra. Although the four observa- 
tions discussed here did not have simultaneous radio spectra, 
our prior studies (Gleissner et al. 2004; Nowak et al. 2005; 
Wilms et al. 2006) allow us to make good estimates of what 
the correlated 15 GHz radio flux likely was. We describe this 
estimated joint radio/X-ray spectrum with the jet model of 
Markoff, Nowak & Wilms (2005) and Maitra et al. (2009). 

It is important to distinguish the jet models as discussed 
in Markoff et al. (2004, 2005) and later works from the 
earlier iteration of the jet model discussed in Markoff et al. 
(2001, 2003). The 2001/2003 work considered a parame- 
ter space where only synchrotron radiation from the jet con- 
tributed to the observed X-rays. Later work, and the version 
of the jet model used here, include X-rays produced not only 
by jet synchrotron radiation, but also by jet synchrotron self- 
Compton processes and Comptonization of seed photons from 
the disk. The synchrotron component is only mildly beamed, 
with peak (3 ~ 2 in the radio emitting portions of the jet. The 
inverse Comptonization component comes from near the jet 
base where /3 ~ 0.4, and thus undergoes very little beaming. 
Thus, although this does mean that the fitted jet parameters are 
dependent upon assumed inclination angle (e.g., the jet input 
power, TVj will depend upon assumed inclination), given the 
low (3 values inherent in the model it is possible to find good 
fit parameters for a wide variety of assumed inclinations. 

A full description of the main jet model parameters can be 
found in the Appendix of Markoff, Nowak & Wilms (2005) 
and in Maitra et al. (2009). The version of the jet model that 
we use is closest to that discussed in Maitra et al. (2009); 
therefore, we repeat part of the model parameter summary 
presented there. We assume a distance of 2.5 kpc and an in- 
clination of 35° for the Cyg X-l system. The jet properties 


are then determined by the input jet power (Ay, expressed as a 
fraction of the Eddington luminosity for a 1OM 0 black hole), 
the electron temperature of the relativistic thermal plasma en- 
tering at the jet base (T e ), the ratio of magnetic to particle en- 
ergy density (the equipartition factor, k ), the physical dimen- 
sions of the jet base (assumed to be cylindrical with radius Ro 
and height ho, the latter fixed to l.SRo), and the location of 
the point on the jet (z acc ) beyond which a significant fraction 
of the leptons are accelerated to a powerlaw energy distribu- 
tion that follows E~ p 9 with p being a fit parameter. The fit 
parameter e acc can be physically interpreted as being related 
to the particle acceleration rate (which is cx e“ c ! c ; see Markoff, 
Falcke & Fender 2001). 

As for the corona-only models, we also include emission 
from a low temperature disk with inner radius R m and peak 
temperature T^iski • (This disk is very similar in temperature 
and normalization to the low temperature disk component in- 
cluded in the Comptonization model fits.) A relativistically 
broadened Fe Ka line (with the same emissivity and inclina- 
tion parameters as for the corona-only models) is coupled to 
the disk component by having its inner radius of emission tied 
to the disk inner radius. A narrow gaussian line is included, as 
well as a reflection component (calculated with the ref lect 
model from XSPEC). Unlike the reflection model internal to 
the eqpair code, the reflect model, which performs a 
convolution on any given input spectrum, does not account 
for relativistic smearing. The Fe edge near 7.1 keV in this 
unsmeared reflected spectrum can produce a sharp feature not 
seen in the data residuals; therefore, we smear this component 
with a unit-normalized, <r = 1 keV gaussian convolution 7 . 

Fit results for this model are presented in Table 6, Fig. 13, 
and Fig. 14. The overall quality of these fits is quite good, 
although not quite as good as for the corona-only models dis- 
cussed above. The jet models leave slightly larger residuals 
in the PC A spectra (near the Fe line region and the lOkeV 
region), as well as in the GSO spectra. The former could be 
related to inadequate modeling of the reflection spectrum (i.e., 
not using a relativistically smeared model). 

Fits to the hardest X-ray spectra (i.e., the cutoff seen in the 
GSO) are governed by the jet SSC and the inverse Comp- 
ton components. This latter component is comprised of a 
magnetized, beamed corona with a high temperature (3- 
5) x 10 10 K, i.e., « 3-5 MeV], and low electron scattering op- 
tical depth (r es S5 0.01). A high temperature, low optical depth 
corona with a limited residence time under those conditions 
simplifies the jet code by allowing one to use a single scatter- 
ing approximation for calculating the Comptonization spec- 
trum. The compactness of such a corona, however, could lead 
to high pair production which would then serve to cool the 
corona and increase its optical depth (see, for example, the 
critique of Malzac, Belmont & Fabian 2009). 

Post-facto estimates of the physical self-consistency of fit- 
ted jet models are discussed by Maitra et al. (2009) who cal- 
culate the pair production and annihilation rates at the base of 
the corona. Fits to the hard state spectra of XTE J1181+105 

7 Kernels for relativistic smearing do exist, e.g., the models of Brenneman 
& Reynolds (2006) and Dauser et al. (2010). The jet model calculations, 
however, are already substantially slower than the Compton corona calcula- 
tions. Coupling them with the computationally expensive relativistic kernel 
would make the fits discussed here prohibitively time consuming to run. Fur- 
thermore, the jet model reflection geometry is likely more complex than that 
of the standard Compton corona: it should have separate contributions for the 
synchrotron, SSC, and Comptonized disk photons (see Markoff & Nowak 
2004). For these reasons we chose a very simple smearing profile. 
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Fig. 14. — Same as Fig. 10, but for the jet model fits. The following model components are shown: the summed model with the reflection component excluded, 
the reflection component by itself, the broad and narrow Fe Ka line components, the jet synchrotron component (which dominates the 2-10keV spectra for 
observations 1,2, and 4), the jet SSC component (which dominates the 1 00-300 keV spectra for observations 1,2, and 4), and the disk component. Neutral and 
ionized absorption are applied to all shown model components. 


yield far higher annihilation rates than production rates, and 
therefore represent a self-consistent solution. Fits to the hard 
state spectra of GX 339-4 yield mixed results. For some 
cases, annihilation and production rates are comparable and 
the resulting corona is marginally self-consistent. In other 
cases, the calculated production rates are an order of magni- 
tude higher than annihilation rates, and the resulting coronae 
are not self-consistent (Maitra et al. 2009). 

We have performed the same calculations as discussed by 
Maitra et al. (2009) for jet models fitted in this work and find 
that all yield pair production rates larger than pair annihilation 
rates. The coronae within these jet models are therefore not 
self-consistent. Observation 3 produces a fit that is closest to 
being self-inconsistent, with pair production and annihilation 
rates being a factor ps 10 apart from one another. Jet fits to the 
other observations are further from self-consistency. 

Such inconsistencies were similarly true for early coronal 
models, e.g., those based upon the work of Haardt & Maraschi 
(1991). A number of such “slab geometry” coronae, in order 
to produce the hardest X-ray spectra observed in sources such 


as Cyg X-l, used coronal temperatures (apart from any con- 
sideration of pair production) that were unachievable in those 
geometries (see the discussion of this issue by Dove, Wilms & 
Begelman 1997). Such inconsistencies were in fact what led 
to the development of more physically self-consistent coronal 
models, e.g., the kotelp model of Stem et al. (1995) and 
Dove et al. (1997) and the eqpair model of Coppi (1999). 

We consider these jet models, however, as they are as of 
yet the only spectral models that make a serious attempt to 
explain the correlated radio spectra. Future iterations of these 
models will incorporate consideration of pair production, as 
well as synchrotron and SSC cooling of the Compton corona 
(Pe’er & Markoff, in prep.). (As shown by Dove, Wilms & 
Begelman 1997 and others, coronal cooling can be a more 
significant consideration than pair production in some situa- 
tions.) 

The derived parameters are similar to those found for fits 
to other BHC (see the fits to XTE J1 118+105 and GX 339-4 
presented by Maitra et al. 2009). The power input at the base 
of the jet, in terms of Eddington luminosity, is comparable to 
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Fig. 15. — Fe Kq line region residuals for observation 4, obtained from our best fits by setting the broad and narrow line normalizations (including the 
normalization of the Fe xxv absorption line) to zero. Left: the high £7^ eqpair fit, the hybrid thermal/non-thermal eqpair fit, and the jet model fit. (Purple 
circles are the RXTE residuals, while the blue histograms are the Suzaku residuals.) 


the inferred bolometric luminosity (i.e., compare the values 
of JVj to the values of Table 2; although different inclination 
assumptions will alter the value of JVj, given the low (3 values 
in the jet, it is always comparable to the inferred system lu- 
minosity). The radius of the base of the jet is comparable to 
the radius of the innermost stable circular orbit (ISCO). The 
largest fit value exceeds the ISCO for a Schwarzschild black 
hole by only a factor of 1 1 . The shock acceleration distance 
along the jet, z acc ranges from » 60-430 GM/c 1 . The former 
value is comparable to that found for XTE J1 1 18+105, while 
the latter is comparable to that found for GX 339-4. As dis- 
cussed above, the fitted coronal electron temperatures at the 
base of the jet range from « (3-5) x 10 10 K. Work in progress 
suggests that a physically self-consistent coronal temperature 
that properly accounts for the dominant effects of synchrotron 
and SSC cooling may truncate the allowed range of coronal 
temperatures to the lower end found here (Pe’er & Markoff, 
in prep.). Future versions of the jet code will explicitly ac- 
count for these effects. 

For all the fits presented in Fig. 14, both synchrotron and 
SSC radiation make substantial contributions to observed X- 
ray spectra. This is an important point to which we shall re- 
turn in §6. 

5.3. Comparison of Implied Relativistic Lines 

In Fig. 15 we show the Fe line region residuals for the two 
Compton corona models and the jet model. These residuals 
are derived by setting the line normalizations (including ab- 
sorption line normalizations) to zero, without refitting. They 
are also different from line residuals that are often used for 
illustrative purposes where only a portion of the continuum 
outside of the line region is fit with a simple powerlaw (e.g., 
Fig. 7). The residuals in Fig. 15 are meant to highlight the 
differences in the derived line parameters for the three main 
continuum models discussed here. 

As shown by the parameters in Tables 4—6, there are in- 
deed differences among the inferred line parameters. For the 
most part, the purely thermal Comptonization model yields 
the strongest line equivalent widths, while the non-thermal 
Comptonization model yields the weakest equivalent widths 
(by up to a factor of two). The thermal Comptonization model 
usually also yields the smallest inner radii for the relativisti- 
cally broadened line, in some cases extending all the way to 
the ISCO for a Schwarzschild black hole. Larger values of 
these radii, however, are typically found with the other two 
models. These differences are evident in the lower energy 
residuals shown in Fig. 7. 


To be clear, a substantially broadened line is required by 
all spectral models that we have explored, and there is a good 
agreement between the PCA and Suzaku -XIS residuals re- 
garding this fact. We also reiterate that for all explored mod- 
els, the equivalent width of the fitted narrow line (typically 
<15 eV) is only a fraction of that for the broad line, and is in 
fact substantially smaller than the Fe line equivalent width we 
have found for any model that we have applied to any RXTE 
observation of Cyg X-l (Wilms et al. 2006). Furthermore, for 
these four observations, we have not found any implied inner 
radius greater than » 40 GM/c 2 . The precise details of this 
broadened line, however, do depend upon the assumed con- 
tinuum model. 

6. DISCUSSION 

In this work we present broad band (0.8-300 keV) fits to 
four separate observations of Cyg X-l that have simultaneous 
Suzaku and RXTE data. The most recent of these observa- 
tions also has simultaneous Chandra-HETG data. Each of 
these observations shows evidence of dipping events likely 
due to dense structures (“clumps”) in the otherwise highly 
ionized wind of the secondary. This is seen in the lightcurves 
(Figs. 1), and the color-color diagrams (Fig. 2, which can 
be modeled by the presence of a dust scattering halo). The 
ionized absorption is very significant in the Suzaku spectra 
(Fig. 5) even though in this work we consider the least ab- 
sorbed periods of the lightcurve. 

• Modeling the ionized line absorption present in Cyg X- 
1 , via the use of parameterized fits to the HETG spectra, 
is crucial for deriving good fits to the soft X-ray spectra 
seen by Suzaku and RXTE. 

We note that even accounting for this ionized absorption, 
the spectral fits presented in this paper yield reduced xt that 
range from 1. 8-2.5. Do such values truly represent good fits 
to these data? Cyg X-l is bright, and these observations are of 
sufficient length that the signal-to-noise values for these spec- 
tra are quite high. The spectra are dominated by systematic 
errors, especially at the soft X-ray energies. We already have 
added 0.5% systematic errors to the PCA spectra, which is a 
reasonable estimate for the internal uncertainty of the PCA. 
There are also relative uncertainties among the detectors. For 
observation 4, if we increase the PCA systematic errors to 1 %, 
add 1% systematic errors to the HEXTE A cluster (i.e., the 
fixed cluster), and add 3% systematic errors to the Suzaku- 
XIS spectra, then the fits presented here would have reduced 
xt ~ 1. Comparing the fits among the individual XIS spectra 
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(representing both different individual detectors and different 
data acquisition modes), we have found that ±3% is a rea- 
sonable estimate of the end-to-end differences among these 
spectra. We hypothesize that the quality of the fits presented 
here is near “optimal” given the current internal and relative 
calibrations of these detectors. 

The observations that occurred on 2007 May 17, i.e., the 
third set of observations, are potentially the most problem- 
atic in turns of cross-calibration issues. For our other sets 
of observations, the RXTE-PCA data act as a “bridge” be- 
tween the Suzaku-XIS and -PIN spectra. Each can have its 
relative normalization anchored by a comparison to the PCA, 
which overlaps the energy coverage of both detectors. This 
is lacking in the third observation, leading to the worry that 
changes in the normalization constant are subsuming, for ex- 
ample, continuum spectra associated with the spectral break 
at « 1 0 keV. For the third observation, we find the ratio of the 
PIN to XIS 1 normalization constant to range from 1 . 1 0-1 .17. 
This is to be compared to the 1 . 1 6-1 . 1 8 value found for fits to 
the first observation, the 1.06-1.08 found for the second ob- 
servation, and the 1 . 1 0-1 . 1 1 found for the fourth observation. 
The expected value 8 for the P1N/XIS 0 comparison is 1.16. 
There is some amount of scatter for the cross-normalization 
values among the different observations; however, the first, 
second, and fourth observations show little scatter in its value 
for different fits to the same data. More scatter is seen for 
the third observation; therefore, additional systematic uncer- 
tainties need to be considered as being present for that set of 
observations. 

Using simple broken powerlaw and exponentially cutoff 
powerlaw fits, we find that these spectra are among the hardest 
seen in the “low hard state” of Cyg X-l over the past decade. 
For all four observations, the spectra are clearly detected out 
to 300 keV with the GSO, and exponential cutoffs are well 
constrained. 

• Although these are among the hardest Cyg X-l spectra 
ever detected, the exponential folding energies vary by 
over a factor 1.5, and range from 160-250keV. 

Historically observed hard state spectra in Cyg X-l have 
shown folding energies that vary over a slightly wider range, 
while hard state BHC as a class show folding energies that 
span a factor of 5. 

As these spectra are among the faintest and hardest for 
Cyg X-l , they make excellent test beds for theoretical models 
that posit, for instance, that the hard state represents a con- 
figuration with an inner disk that has evaporated into a quasi- 
spherical corona. In such a scenario, we might expect these 
spectra to show the most “extreme” recession of the inner 
disk, although we have noted that the bolometric luminosi- 
ties represented by these spectra only span a factor of two. At 
a few percent of the Eddington luminosity, they are not far 
below the expected soft-to-hard state transition. Numerous 
transient BHC sources show much fainter hard states as they 
fade into quiescence. 

We have presented a number of different spectral models, 
all of which describe the 0.8-300 keV spectra well. Some of 
these models describe the X-ray spectra primarily with Comp- 
tonization components (whether due to a thermal or hybrid 
thermal/non-thermal corona), while the jet model is domi- 

8 http : //heasarc. gsfc.nasa . gov/docs/suzaku/analysis/ 
watchout . html . Since we did not fit XIS 0 data for the fourth observa- 
tion, here we compare to XIS 1. 


nated by synchrotron and SSC emission from the jet. All of 
these models have a number of features in common. 

• All models require a soft excess that here we describe 
with a disk component with low (AJdisk ~ 200 eV) peak 
temperature. The implied inner radii of these disks 
range from 2^40 GM/c 2 . 

• All models require a relativistically broadened line 
component. The inner emission radius of this broad- 
ened line never exceeds « 40 GM/c 2 , but for some 
models is as low as 6 GM/c 2 . 

• All models require a reflection component. The typical 
values for the reflection fraction are D/27T ~ 0.2-0.3. 

• All models imply that the spectral hardening at « 
10 keV is not solely due to reflection. 

This latter point is very important, and broadly agrees with 
similar conclusions drawn by Frontera et al. (2001), Ibragi- 
mov et al. (2005), and Makishima et al. (2008). The presence 
of the broad Fe line and some of the spectral curvature in the 
20-300 keV band is a clear indication of the presence of re- 
flection. However, this reflection spectrum is not sitting on 
top of a simple “disk+powerlaw” spectrum. There is ad- 
ditional continuum complexity separate from reflection that 
contributes to this perceived break. A high seed photon tem- 
perature in the thermal corona model yields a soft excess in 
the 2-10 keV band and thus contributes to the measured break 
in that scenario. (See the discussion in Wilms et al. 2006.) As 
has been discussed by Ibragimov et al. (2005), a non-thermal 
electron population in the corona can lead to a soft excess 
in the 2-10keV band that helps contribute to what otherwise 
would be modeled as a reflection break at lOkeV. The two 
corona model discussed by Makishima et al. (2008), a ver- 
sion of which is shown in Fig. 12, rather explicitly replaces 
part of the reflection component with a broad band contin- 
uum model. Finally, in the jet paradigm the spectral break at 
lOkeV is partly attributable to the transition from dominance 
by synchrotron emission to SSC emission in the continuum. 

There are plausible physical scenarios for each of the dis- 
cussed spectral models, with some hypothesized geometries 
being shown in Fig. 9. The hybrid thermal/non-thermal coro- 
nal model is the closest to the concept of the quasi-spherical 
inner corona with outer geometrically thin disk (e.g., Eard- 
ley, Lightman & Shaipro 1975; Shapiro, Lightman & Eardley 
1976; Ichimaru 1977; Dove et al. 1997, etc.). We have only 
been able to find such solutions, however, when invoking a 
partly non-thermal electron population in the corona. 

The purely thermal coronal model contains two disk com- 
ponents that are reminiscent of the physical description given 
by Mayer & Pringle (2007). These authors describe a situa- 
tion where an outer, geometrically thin, cool disk surrounds 
an inner, geometrically thick, hot corona. In the very inner 
radii of this corona, however, thermal conduction leads to it 
condensing into a geometrically thin and optically thick disk 
(see Fig. 9). Such a component could supply the high tem- 
perature seed photons in our thermal corona solutions, while 
the lower temperature outer disk could provide the bulk of the 
reflection features. The jet model has a natural physical in- 
terpretation in that the usually observed optically thick radio 
spectrum observed in the hard state is clear indication of the 
presence of a jet. The question that remains is the contribution 
of this component to the X-ray band. 
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The fact that the continuum is more complex than a simple 
“disk+powerlaw”, yet there are multiple, physically mo- 
tivated models that yield comparably good spectral fits, leads 
to the final point. 

• Although a relativistically broadened line is required in 
all of our spectral models, the parameters of this line 
are dependent upon the presumed continuum model. 

Coupled with this dependence upon assumed continuum spec- 
trum is an implicit dependence upon ionized absorption, for 
which we have detailed Ch an dra-HETG measurements for 
only the fourth observation 9 . Again, we have not found an 
inner radius for this line that exceeds « 40 GM / c 2 . Given the 
variations of this line with presumed continuum, however, we 
are as of yet unable to use this line for more refined diagnos- 
tics such as estimates of black hole spin. 

Although we are unable as of yet to draw firm conclusions 
as to the best geometry and physical mechanisms to describe 
the hard state spectra of Cyg X-l, these new joint Suzaku - 
RXTE data provide a stunning contrast to our prior results 
using solely RXTE data (Wilms et al. 2006). For the RXTE 
data alone we were able to describe the 3-125 keV spectra 
with a variety of physically motivated Comptonization mod- 
els and to describe correlations among the fit parameters. On 
the other hand, the simple exponentially cutoff, broken pow- 
erlaw models with a single, broad gaussian line described the 

9 Since performing these observations, we have carried out a Chandra - 
HETG observation of orbital phase 0.5 (PI: Nowak), and have an approved 
observation of orbital phase 0.25 (PI: Hanke). Coupled with archival Chan- 
dra-HETG observations of orbital phases near 0.75, we hope to develop a bet- 
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the sophistication and physical self-consistency of the spec- 
tral models will continue to improve. Together, they offer the 
promise of obtaining a better understanding of the physical 
processes occurring in these BHC systems. 
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Parameters for thermal Comptonization spectral fits 
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Note. — See text for a description of the model. Spectra are normalized to PCA data. Disk normalizations are (/?in/D) 2 cos0, where R\n is the inner disk radius in units of km, D is the source distance in units 
of lOkpc, and 0 is the inclination. Line normalizations are integrated photon flux in the lines ( 7 /cm 2 /s). Error bars are 90% confidence level for one interesting parameter. 
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Parameters for non-thermal Comptonization spectral fits 


Date 

(yyyy-mm-dd) 

Nr 

( 10 22 cm- 2 ) 

^eqp 

4 /is 

4t / 4 

Tp 

Q/27T 

$ 

(47Tfion/n) 

/?in 

(GM/c 2 ) 

/-/diskl 

(keV) 

diskline 

( 10 - 3 ) 

-<4 line 
( 10 ^) 

Qjnes 

X^/DoF 

2006- 10-30 

2007- 04-30 

2007- 05-17 

2008- 04-19 

0.96/gl 

°' 86 S? 

i- 08 X 

9 99Q+0.027 
-Q-QS2 
1 S8+0003 

1 33R+8:87? 

16 -04-8 : ?3 

13-29 

12 - 98 5? 

15-28t| 

0-01 

°- 36 5S 

o- 78 3:i2 

1 tq+0.03 
-Q-Q? 

0 4o+# 0 1 

0 01 + ^ 
-QQI 

1 41+0.01 

- 0.02 

0 21 l+O-OOi 
-0.008 
0 305 + ^^ 

H 9 S 7+0 0 i 8 
- 0.002 

499 / 33 

97+i° 

0+ 18 

'89! 43 

6 - 7 l! 

13.3/0;} 

6 3 +20 
-0,3 
09 0 + 1.8 
JZ “ J -0.2 

O 1 Q 9 + 0.002 

-Q-QQi 

O 1 S 8 + 0-001 

U '° -Q-QQi 
0 175+0003 
• -Loon? 

0 150 4 *™ 

u.ijUq 001 

6 . 6 ^| 
3 O 4 ^ 6 

41^ 
H - 1.0 

11 1 +31 

- 1-2 
2 6 +1 - 8 
94+11 

-4.4 

a32 H 

O-SS/g 

n S1+0O7 

1 00 + ° ° 7 
1 - uu _0.07 

2254/977 

1919/844 

810/485 

1535/845 


Note. — See text for a description of the model. Spectra are normalized to PCA data. Disk normalizations are (R^/D) 2 cos0, where R ^ is the inner disk radius in units of km, D is the source distance in units 
of lOkpc, and 0 is the inclination. Line normalizations are integrated photon flux in the lines ( 7 /cm 2 /s). Error bars are 90% confidence level for one interesting parameter. 
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Parameters for jet-dominated spectral fits 
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Note. — See text for a description of the model. Spectra are normalized to PCA data. Line normalizations are integrated photon flux in the lines ( 7 /cm 2 /s). Error bars are 90% confidence level for one 
interesting parameter. 
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APPENDIX 

Suzaku ATTITUDE CORRECTION AND PILE-UP ESTIMATION 

We have created two tools to aid our Suzaku data analysis: aeattcor . si and pile_estimate . si. The tools, descrip- 
tions of their use, and example results can be found at: http : //space .mit . edu/ASC/software/suzaku. The former 
tool further corrects the Suzaku attitude solution, whereas the latter estimates the degree of pileup in a given subset of the obser- 
vation. Both are scripts written in S-Lang but are designed to be run on a terminal command line using ISIS, the Interactive 
Spectral Analysis System (Houck & Denicola 2000), as a driver. For all intents and purposes, the scripts behave similar to the 
typical Unix command line tools found in the HEASOFT package. 

Thermal flexing of the Suzaku spacecraft leads to a slow wobbling of the optical axis, and hence blurring of the image. 
Current Suzaku tools partially correct this effect by adjusting the spacecraft attitude file based upon details of the spacecraft 
orbit, temperature, etc. (Uchiyama et al. 2008). aeattcor . si further improves this correction (Fig. 16) for a bright source by 
using its time-dependent detected image to create a new attitude file. Specifically, the tool attempts to shift the time-dependent 
mean detector image position to a new specified, fixed sky position. It presumes that there are no intrinsic variations in the time- 
dependent image position. (A mean image position may be time varying, for example, if it is comprised of two or more variable 
sources of comparable flux. Additionally, a highly variable and piled up source where an “image crater” appears and disappears 
over the course of the observation may lead to a variable mean image position.) The bin time over which the tool searches for 
and shifts the image peak is a user selectable parameter; however, we have used the default of 100 seconds. 

We applied the tool to detector images that first underwent the standard attitude correction process (Uchiyama et al. 2008). We 
chose circular regions with radius 2' that were visually centered on the source image. The tool then created a new attitude 
correction file, which we applied with the xiscoord tool. All of our attitude correction files were created using solely the XIS 
1 images. 

Pileup occurs in CCD detectors when two or more photons fall on the same or neighboring pixels during the same 
readout frame, and therefore are read as a single higher energy photon or are discarded as a bad event (Davis 2001). 
pile_estimate . si was run after aeattcor . si, as uncorrected blurring leads to an underestimate of the degree of pileup. 

The tool first creates a lightcurve in three energy bands, using both rate and counts. Here we used the default values of 0.5- 
1 .5 keV, 1 .5-3 keV, and 3-9 keV for the energy bands and 32 seconds for the time bins. These data are then passed to the vwhere 
(Noble 2005) filtering tool. The vwhere tool was used to create the color-intensity diagrams shown in Fig. 2), and it was used 
to select times of approximately uniform rates and colors (see §3). These time selections were then written to a filter file used for 
subsequent data extraction using the xselect tool. 

A box-car smoothed (3x3 pixel bins) image is displayed using the ds 9 tool, with the image rescaled to the estimated pile up 
fraction. Pile up fraction here is defined to be the ratio of events lost via grade or energy migration to the events expected in the 
absence of pileup. Furthermore, the pile up fraction is based upon the mean counts per 3x3 pixel region per readout frame. The 
image displays discrete steps that represent the minimum pileup fraction in the displayed region, with the exception of the highest 
displayed value. This value corresponds to the image maximum, and is shown over regions at this maximum down to half way 
towards the next lowest displayed value. 

Using this image, we chose to exclude the central regions with pile up fractions > 10% (Fig. 16). Typically, this was an « 20" 
radius region, comprising « 1/3 of the detected events. The pile_estimate . si tool then gives an estimate of the mean level 
of the remaining pile up fraction in the image, which for these data was typically < 4%. 
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Fig. 16. — Left: An example of Suzaku attitude correction. The left half of the image shows a Cyg X-l image using the standard attitude correction. The right 
half of the image shows the improvement with aeattcorr.sl. Right: An image of Cyg X-l where discrete colors correspond to the pileup in that region. The 
outer white circle denotes the outer boundary of an annular extraction region, and the inner white circle denotes the inner boundary. Pileup fractions within this 
excluded region are as high as 35%. The average effective residual pileup level is < 4%. 





